leaflet的小搭档leaflet.minicharts来了，从此动态地图又多了一些乐趣~~~
本文内容取材自leaflet.minicharts包官方主页的案例介绍，本篇案例虽然是关于leaflet在线地图的辅助包，但是该包的出现对于leaflet生态系统来说，确是有着划时代的意义。
该包大大扩充了leaflet包所能呈现的图表形式，打破了散点图、路径图、热力图三类图表对于传统地图数据呈现形式的垄断地位。
该包为此提供了气泡饼图、玫瑰图、mini柱形图以及时间维度控制器等多种数据可视化元素。
以下便是作者对该包的简要介绍及案例演示：
一直以来，借助于Rstudio团队开发的交互式地图工具包——leaflet，在R语言中创建交互式地图已非难事。但是该包仅仅提供了很少的集中图表呈现类型，因而在一张单独的地图上呈现的数据信息量被大大限制了。
倘若你要在地图上呈现一些点数据，你只能需要使用两个变量定位出这些点在地图上的位置，然后更改点的半径和点的填充颜色来呈现数据。
然后leaflet.minicharts包的出现大大改变了这一格局。你可以在利用其提供的两个附加函数，在leaflet包的交互地图上增加更多的mini图表。
这些mini图表可以呈现比之前多得多的变量数据（理论上不受变量数量的限制），目前为止，仅支持三种图表类型：
barcharts (the default), pie charts and polar area charts.
该包内提供了一个法国自2010~2017年包含电力产品、电力消费及结构变动趋势的数据集。
library("ggplot2")
library("leaflet.minicharts")
library(geojsonio)
library(rgdal)
data("eco2mix")
head(eco2mix)
geojson <-readOGR("D:/R/mapdata/State/France.json","OGRGeoJSON",stringsAsFactors=FALSE)
#导入json格式的法国地图数据素材
library(dplyr)
prod2016 <- eco2mix %>%
mutate(
renewable = bioenergy + solar + wind + hydraulic,
non_renewable = total - bioenergy - solar - wind - hydraulic
) %>%
filter(grepl("2016",month) & area!= "France") %>%
select(-month) %>%
group_by(area,lat,lng) %>%
summarise_all(sum) %>%
ungroup()
#抽取了法国2016年的能源消费结构数据信息：
library(leaflet)
library(leafletCN)
basemap <- leaflet(geojson,width ="100%",height="400px") %>% amap() %>%
addPolygons(weight=1,color="#444444",opacity=1,fillColor ="white",fillOpacity =0.7,smoothFactor=0.5)
#将地图背景保存为临时对象：
renewable2016 <- prod2016 %>% select(hydraulic,solar,wind)
colors <- c("#3093e5","#fcba50","#a0d9e8")
basemap %>%
addMinicharts(
prod2016$lng, prod2016$lat,
chartdata = renewable2016,
colorPalette = colors,
width = 45, height = 45
)
#默认图表形式是柱形图
#气泡饼图
colors <- c("#4fc13c", "#cccccc")
basemap %>%
addMinicharts(
prod2016$lng, prod2016$lat,
type = "pie",
chartdata = prod2016[, c("renewable", "non_renewable")],
colorPalette = colors,
width = 60 * sqrt(prod2016$total)/sqrt(max(prod2016$total)),transitionTime=0)
#玫瑰图1：
basemap %>%
addMinicharts(
prod2016$lng, prod2016$lat,
type = "polar-radius",
chartdata = prod2016[, c("renewable", "non_renewable")],
colorPalette = colors,
width = 60 * sqrt(prod2016$total)/sqrt(max(prod2016$total)),transitionTime=0)
#极坐标（玫瑰图）
#玫瑰图2：
basemap %>%
addMinicharts(
prod2016$lng, prod2016$lat,
type = "polar-area",
chartdata = prod2016[, c("renewable", "non_renewable")],
colorPalette = colors,
width = 60 * sqrt(prod2016$total)/sqrt(max(prod2016$total)),transitionTime=0)
添加时间趋势：
#柱形图（带时间趋势）
prodRegions <- eco2mix %>% filter(area != "France")
basemap %>%
addMinicharts(
prodRegions$lng, prodRegions$lat,
chartdata = prodRegions[, c("hydraulic", "solar", "wind")],
time = prodRegions$month,
colorPalette = colors,
width = 45, height = 45
)
#饼图（带时间趋势）
prodRegions <- eco2mix %>% filter(area != "France")
basemap %>%
addMinicharts(
prodRegions$lng, prodRegions$lat,
chartdata = prodRegions[, c("hydraulic", "solar", "wind")],
time = prodRegions$month,
colorPalette = colors,
type="pie",
width = 45, height = 45
)
#玫瑰图（带时间趋势）
prodRegions <- eco2mix %>% filter(area != "France")
basemap %>%
addMinicharts(
prodRegions$lng, prodRegions$lat,
chartdata = prodRegions[, c("hydraulic", "solar", "wind")],
time = prodRegions$month,
colorPalette = colors,
type="pie",
width = 45, height = 45
)